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ABSTRACT 



A DIRECT METHOD FOR THE MEASUREMENT 
OF CONTROL ROD WORTH IN A NUCLEAR REACTOR 

By 

William Arthur Emslie 



It is often necessary to perform control rod calibration 
in a nuclear reactor. Such calibration requires the measurement 
of the reactivity produced when discrete sections of a control rod 
are withdrawn from a critical reactor pile. From these data a 
total reactivity worth curve can be complied for the length of the 
rod. The subject of this thesis is the development of a means by 
which reactivity worth can be measured. 

A solution to this measurement problem has been developed 
using the current output of a compensated ionization chamber 
(adjacent to the reactor) as a signal source. When a section of 
control rod is withdrawn from a critical reactor, the current from 
the ion chamber is 



i 

s 



(t> 



= V 



t/9 



P 



where t is time, 0^ is the reactor stable period and 1^ is 
the initial ion chamber current. 

The stable period is proportional to the amount of reactivity 



produced since 



William Arthur Emslie 



0 o e o 
e p - — a T 

ex 

is the mean neutron lifetime, k is the excess multiplication 
0 ex 

factor and p is the reactivity. This approximate relationship 
between 0^ and p makes it possible to obtain the reactivity 
from the output current of the ion chamber since 



yt) = I Q e 



t/e 



pt/e 



- I o e 



Actual measurement of p is accomplished with a circuit 
composed of a logarithmic amplifier followed by a differentiator 
and filter. The transfer function of this circuit is 



(-R^s) x 

V out (s) = V l (s) (R^s+1) (R 2 C 2 s+ 1) (R f C f s+l) 



where V^(s) is the frequency domain expression for V^(t) (the 
output voltage of the logarithmic amplifier) and is equal to 



V^s) = 



~ A 1 ° S 10 ( - I 0^ I k' ) A log 10 e 

S 2 

9 p s 



In the time domain, the logarithm of the ion chamber current 



produces a linear ramp voltage v^(t) where 



t/0 



v l (t) = -A log 10 



i <t> 

s 



= -A log 



V 



10 



AT T 0 At 1 

- A log io r - t log io e 
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This ramp voltage is then differentiated and filtered to 
obtain a d-c voltage proportional to the inverse of p. The re- 
sulting voltage (V^) is 
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V 2 ‘ - K k <V e » * - K i <-* log 10 if - f 



l°g 10 e > 
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V 2 may be calibrated on a d-c voltmeter to indicate the 
corresponding reactivity. 

A circuit to synthesize the transfer function was designed 
and constructed. Results of tests made on circuit performance show 
the total accuracy in the measurement of p to be within +0.5 
percent of the meter scale. 



A DIRECT METHOD FOR THE MEASUREMENT 
OF CONTROL ROD WORTH IN A NUCLEAR REACTOR 

By 

William Arthur Emslie 



A THESIS 
Submitted to 

Michigan State University 
in partial fulfillment of the requirements 
for the degree of 

MASTER OF SCIENCE 

Department of Electrical Engineering and Systems Science 



1972 



Tr\* s 
F a* ^ f 






TO MY PARENTS 



ii 



ACKNOWLEDGMENTS 



The author wishes to express his sincere thanks to his 
major professor. Dr. G.L. Park, for his guidance and support during 
the course of the research. Sincere appreciation is also expressed 
to the other members of the guidance committee, Dr. B. Wilkinson 
for his advice on nuclear reactor theory and Dr. L. Giacoletto for 
his guidance on the design of the reactivity measuring circuit. 

Special thanks is extended to Mr, E. Brockbank who made 
possible extensive testing and calibration of the reactivity meter. 



iii 



TABLE OF CONTENTS 



Page 

LIST OF TABLES v 

LIST OF FIGURES vi 

LIST OF SYMBOLS viii 

Chapter 

I. INTRODUCTION 1 

II. REACTOR TRANSIENT BEHAVIOR 5 

2.1 Reactivity, Reactor Period and Theoretical 

Relationships 5 

2.2 Relationship Between Reactivity and the 

Reactor Stable Period 7 

2.3 Effects of Prompt and Delayed Neutrons on 

Reactor Transient Behavior 7 

2.4 The Compensated Ion Chamber H 

III. DESIGN OF THE REACTIVITY MEASURING DEVICE I 5 

3.1 Formulation of the Problem 15 

3.2 Possible Methods of Solution 16 

3.3 The Circuit and Theory of Operation 20 

3.4 Calibration 23 

IV. CONSTRUCTION OF THE REACTIVITY METER .... 27 

4.1 Circuit Tests, Performance and 

Specifications 27 

V. CONCLUSIONS 30 

BIBLIOGRAPHY .... 32 

APPENDIX A DESIGN OF THE LOGARITHMIC AMPLIFIER 33 

APPENDIX B DETERMINATION OF THE RAMP FREQUENCY COMPONENTS 3 9 

APPENDIX C DESIGN OF THE DIFFERENTIATOR 42 

APPENDIX D DESIGN OF THE FILTER 47 

APPENDIX E SPECIFICATIONS OF THE COMPENSATED ION CHAMBER . 49 

iv 



LIST OF TABLES 



Table Page 

4,1 Specifications of the reactivity meter 29 

A.l Electrical specifications of the Burr-Brown 

modular logarithmic amplifier model 4116 (current 
input only) 34 

E.l Specifications of the Westinghouse type WL-8501 

compensated ion chamber 49 



v 



LIST OF FIGURES 



Figure Page 

1.1 Calibration curve for the shim rod in the 
Michigan State University TRIGA Mark I 

nuclear reactor 3 

2.1 Relationship between reactivity and reactor 

stable period for the Michigan State University 
TRIGA Mark I nuclear reactor (Courtesy of Gulf 
General Atomic) 8 

2.2 Combined effect of prompt and delayed neutrons 

in a supercritical reactor 10 

2.3 Simplified diagram of a compensated ion chamber . 12 

2.4 Graph of ion chamber output current as a function 
of reactor output power level for the Michigan 

State University TRIGA Mark I nuclear reactor ... 14 

3.1 Block diagram of a circuit which provides a d-c 
voltage proportional to the inverse of the 

reactor stable period 16 

3.2 Circuit for obtaining a d-c voltage proportional 

to 0 17 

P 

3.3 Circuit for obtaining a discrete voltage 

proportional to the slope of i (t) 18 

s 

3.4 Circuit diagram of the reactivity measuring 

instrument 21 

3.5 Apparatus used to calibrate the reactivity meter 24 

3.6 Graph showing the relationship between reactivity 
and reactor stable period for the Michigan State 
University TRIGA Mark I nuclear reactor 

(Courtesy of Gulf General Atomic) 26 

4.1 Photographs of the reactivity meter 28 

A.l Burr-Brown modular logarithmic amplifier model 

4116 with external connections 35 

vi 



Figure Page 

A. 2 Design curves for the Burr-Brown model 4116 

log amplifier, (a) Reference current vs . R , 

(b) Voltage constant vs. R^ ¥ 37 

C.l Transfer function plot for an ideal 

differentiator, (a) Imaginary part, (b) Real 

part 42 

C.2 Transfer function plot for a compensated 

differentiator, (a) Imaginary part, (b) Real 

part 43 

C. 3 Equivalent circuit of a compensated 

differentiator 44 

D. l Output filter for the reactivity meter 47 



vii 



LIST OF SYMBOLS 



A Voltage constant for log amplifier 

A (u,) Open loop voltage gain 

C Capacitance 

Iq Initial ion chamber output current 

1^ Reference current for log amplifier 

I Ion chamber output current 

k „ Reactor effective multiplication factor 

eft 

k Reactor excess multiplication factor 

ex 

n Neutron 

P Reactor power level 

Pq Initial reactor power level 

R Resistance 

t T ime 

V Voltage 

Z Impedance 

(3 Beta particle 

(3^ Fraction of delayed fission neutrons in the ith group 

y Gamma photon 

Qq Mean neutron lifetime 

Reactor power level doubling time 

0 Reactor stable period 

P 

9 t Minimum ramp "on" time 

viii 



Decay constant of ith group 
Excess reactivity (reactivity) 
Reactor neutron flux 
Initial reactor neutron flux 



Radian frequency 



CHAPTER I 



INTRODUCTION 



The necessity to perform calibration of the control rods 
in a nuclear reactor requires a means by which to measure the 
reactivity produced by withdrawal of discrete sections of a con- 
trol rod from a critical reactor pile. 

Use of an ionization chamber adjacent to the reactor pro- 
vides an exponential current (i^(t)) from which the reactivity 
may be derived. When a section of control rod is withdrawn from 

a critical reactor i (t) becomes 

s ' 



where 
is the 



i 

s 



(t) 



= I o e 




Iq is the initial ion chamber current, t is 
reactor stable period.'*’ 

The reactivity may be obtained from equation 



( 1 . 1 ) 



t ime and 0 

P 



(1.1) since 



9 



P 




where 



Qq is the mean neutron lifetime. 



(1-2) 



The reactor stable period used in this expression is the time 
required for the ion chamber current (i (t)) to increase by a 
factor of e. 

2 

This approximation is explained and justified in Section 2.1. 



1 



2 



The relationship between p and i (t) expressed in 

s 

Equations (1.1) and (1.2) makes it possible to obtain p from 

i (t) . This may be accomplished by means of a circuit which first 
s 

takes the logarithm of i^(t) and then differentiates and filters 

the result to yield a d-c voltage proportional to p. 

An example of a control rod calibration curve is illustrated 

in Figure 1.1. The curve is compiled by withdrawing the rod a 

short distance and measuring the corresponding reactivity produced. 

This procedure is repeated until the entire length of the rod has 
3 

been withdrawn. 

Chapter II of this study details the theory underlying 
reactivity and the ionization chamber. The effects of prompt and 
delayed neutrons are also discussed. 

In Chapter III the reactivity measurement problem is 
formulated with a set of possible solutions. A particular solution 
method from these alternatives is then presented. Analysis of the 
design, mathematical properties, and calibration of the circuit follows. 
Further design analysis is presented in the Appendices. 



The units of reactivity as expressed in Figure 1.1 are in dollars 
and cents where 1 dollar = 100 cents and is defined as 

k 

1 dollar = -S* ~ £ (1.3) 

0. p. 

*1 

where k is the reactor excess multiplication factor and ft. is 
ex K M i 

a normalizing constant equal to the total fraction of delayed neutrons 

and varies with different reactors. For the Michigan State University 

TRIGA Mark 1 reactor 0.0073 and therefore one dollar's worth 

of reactivity is that amount required to produce an excess 

multiplication factor of 0.0073. 



React ivity (dollars) 



3 




Rod Position (arbitrary units) 

Figure 1.1 Calibration curve for the shim rod in the Michigan 
State University TRIGA Mark I nuclear reactor 
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Chapter IV describes a prototype instrument which was built 
from the design in the previous chapter. Test results and circuit 
specifications are also given. 

Conclusions are made in Chapter V. 



CHAPTER II 



REACTOR TRANSIENT BEHAVIOR 



2.1 Reactivity, Reactor Period and Theoretical Equations 



The excess reactivity (hereafter referred to as reactivity) 
in a nuclear reactor is defined as the ratio of excess to effective 
multiplication factors and is given by 



P = 



ex 



ef f 



( 2 . 1 ) 



In the above formula the effective multiplication factor. 



is equal to unity for a critical reactor. The excess multi- 



plication factor is 



k = k 
ex ef f 



( 2 . 2 ) 



and is equal to zero for a critical reactor. 

Substitution of Equation (2.2) into Equation (2.1) yields 



P 




(2.3) 



Since the values of k __ used in this analysis are to be no 

eft 



The content of this section, except for that material noted 
otherwise, is taken from M.M. El-Wakil, Nuclear Power Engineering 
(New York: McGraw-Hill Book Co., 1962), pp. 128-130. 



3 
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greater than 1.03^, Equation (2.3) may be approximated as 



P = k eff 



1 - k 



ex 



(2.4) 



Using the above approximation for reactivity it is possible 
to relate p to the instantaneous neutron flux within the reactor. 

The equation describing the flux in a supercritical reactor 
is 

t/e 



$ = $ Q e 



(2.5) 



The reactor stable period (0 ) may be written as 

e ' 

P k 
^ ex 

and substitution of Equation (2.4) for k relates n to Q 

ex r p 

where 



(2.6) 



e ^ 

P P 



( 1 . 2 ) 



This approximate expression for 0^ may be substituted 
into Equation (2.5) to form a relationship between reactor flux 
and reactivity such that 

pe/e 0 

S a § 0 e “ (2.7) 



Reference to the rod calibration curve in Figure 1.1 shows the 
maximum reactivity to be approximately four dollars. This value 
corresponds to an effective multiplication factor of 1.0292 as 
defined in the Nuclear Reactor Operations and Training Manuel , 
Reactor Theory (Michigan State University, 1969), pp. 4-5. 



7 



2.2 Relationship Between Reactivity and the Reactor Stable Period 

As Equation (1.2) illustrates, the approximate relation- 
ship between the reactor stable period and the reactivity is weighted 
by the mean neutron lifetime. This equation is characterized by a 
nonlinear property of 0^ which is not constant but rather varies 
slightly with p. Examination of the exact expression for reactivity 
reveals the source of this nonlinearity. p may be written exactly as 



P = 



9 k 

p eff 



2-^ 
v 1 + \, 



\ . 9 
i P 



( 2 . 8 ) 



where 0^ is the fraction of delayed fission neutrons in the ith 

group and \ is the respective decay constant.^ 

For large 0 , Equation (2.8) may be approximated as^ 

P 



0 



, !_ „ h _ 9 

p ^ Zj 

9 p i h % 



(2.9) 



This approximation, however, is only valid for 0 > 100 and as 

P 

0p decreases in value both the first term on the right in Equa- 
tion (2.8) and the 1 in the denominator of the summation no longer 
remain negligible. The net effect is a nonlinear relationship 
between 0^ and p as illustrated in Figure 2.1. 



2.3 Effects of Prompt and Delayed Neutrons on Reactor Transient 
Behavior 



The fission process in the Michigan State University TRIGA 
Mark I nuclear reactor is described by the following nuclear 

Samuel Glasstone, Nuclear Reactor Engineering (New York: 

D. Van Nostrand Co. Inc., 1963), pp . 244-245. 

^ lb id . , pp . 244-245. 



p, Reactivity (cents) 
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0 , Period 
P 



(seconds) 



Figure 2.1 Relationship between reactivity and reactor stable 

period for the Michigan State University TRIGA Mark I 
nuclear reactor (Courtesy of Gulf General Atomic) 
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equations : 



..235 . 1 m 236 j. „ lm 

92 U + 0 n "* 92 U + ener Sy (2.10) 



or 



235 1 1 

92 U + Q n ^ fips + energy + y + c^n (2.11) 

Equation (2.10) occurs about twenty percent of the time. The 

fission products (fips) in Equation (2.11) have a variety of atomic 

235 

weights but will generally weigh about half as much as the U 
Heat is generated by the energy given off in both reactions and the 
energy is about 200 Mev per fission. On the average between two and 
three (2 ^ c £ 3) neutrons are given off per reaction in Equation 
(2.11). This equation is critical to reactor operation since the 

g 

neturons produced must sustain the fission chain reaction. 

Neutrons emitted in Equation (2.11) are of two types: 

Prompt neutrons and delayed neutrons. Prompt neutrons make up 
about 99.27 percent of all neutrons emitted in the TRIGA reactor. 

The total mean lifetime of these prompt neutrons is less than one 
millisecond. ^ 

Delayed neutrons are the result of the decay of certain 
fission products in Equation (2.11). An example of a delayed 
neutron emission is given by the nuclear equation 



53 



x 137 p x 137 



n 



22 sec 



54 e 



ins tantane ous 



v 136 

54 e 



( 2 . 12 ) 



9 



Training Manuel, Reactor Theory , p. 1. 
Ibid . , p. 4. 
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where 



I 



137 



is one of the fission products. The 22 second delay 



53 



in the neutron emission gives it a half life of appearance of 22 

, 10 

seconds . 

Even though the fraction of delayed neutrons is small, their 
long apparent lifetimes make the average lifetime of all neutrons 
emitted much longer than the lifetime of prompt neutrons alone. 

This effect is vital since a reactor period caused by prompt neutrons 
alone would make control of the fission process impossible. 



neutron flux can be described as follows: In a critical reactor 

withdrawal of a control rod some discrete distance will produce a 
transient period due to prompt neutrons and a stable period due to 
the combined effect of prompt and delayed neutrons. The effect of 
the transient period will result in a "bulge 11 in the neutron flux. 
This effect is shown graphically in Figure 2.2."^ 



time (seconds) 

Figure 2.2 Combined effect of prompt and delayed neutrons in 
a supercritical reactor 



The effect of prompt and delayed neutrons on the change in 




p positive 



p = 0 



10 



El-Wakil, Nuclear Power Engineering , p. 131. 
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Ibid . o pp. 133-134. 
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Since p is proportional to the slope of the curves in 
Figure 2.1 it is necessary to let the transient period die out 
prior to measurement of p. 

2.4 The Compensated Ion Chamber 

It has been stated in Equation (2.7) that the reactivity 
produced when a discrete portion of control rod is withdrawn from 
a critical reactor pile is exponentially proportional to the change 
in neutron flux. In order to measure p it is therefore necessary 
to monitor this flux. This is accomplished in the TRIGA Mark I 
reactor by means of the output current of a compensated ion chamber 
located near the reactor core. 

The output of the ion chamber is a very low current which 
is directly proportional to the neutron flux in Equation (2.5) 
such that 

t/e 

= V P 

where 0 may be substituted by the approximate expression 




as previously stated. 

Figure 2.3 is a simplified diagram of the compensated ion 
chamber. The device is actually composed of two chambers where 
R g is the equivalent output resistance. One chamber is lined 
with boron enriched with and is biased at + 580 volts. From 
this sub-chamber a current proportional to neutron flux and gamma 
flux is obtained. The other sub-chamber has a variable bias from 



12 



0 to -37 volts and is sensitive to gamma flux only. If the gamma 

sensitive current is subtracted from the gamma, neutron current 

of the boron coated sub -chamber, a current proportional to neutron 

12 

flux alone is obtained. 



580 V 



Boron 

Coating 



< 



R 



I 



Equivalent Output 
Res is tance 



-37 £ V <; 0 V — , 

c + 



Figure 2.3 Simplified diagram of a compensated ion chamber 

2 2 

At low neutron flux levels ($ < 2.5 X 10 n/cm /sec), the 

current contribution from the gamma sensitive chamber is nearly 

equal to that of the gamma, neutron sensitive chamber. It is 

therefore important to set the compensating chamber voltage (V^) 

to the proper value in order to achieve a linear response due to 

neutron flux alone. This leaves room for error in the linearity 

of i g (t) if is not properly adjusted. Fortunately, as the 

neutron flux level increases the current contribution to i (t) 

s 



12 



Training Manuel, Instrumentation , p. 5. 
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due to gamma flux alone becomes negligible and i (t) approaches 

s 

a linear dependance on neutron flux alone. 

Experimental results reveal that improper adjustment of 

the compensating voltage at low reactor flux levels (corresponding 

to reactor output power levels below ten watts) causes the output 

current of the ion chamber (i^(t)) to become nonlinear. Figure 

2.4 shows the maximum nonlinear effects at low power levels when 

the compensating voltage is adjusted to zero volts. 

Since a linear relationship between ion chamber current 

(i g (t)) and neutron flux (also reactor output power level) is 

necessary, it is best to rely only upon output currents which 

correspond to reactor power levels greater than ten watts. 

Specifications of the ion chamber are given in Appendix E. 

13 

Since the output resistance is greater than 10 ohms, the chamber 
may be considered as a current source. 



Ion Chamber Output Current (nA) 
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Figure 2.4 Graph of ion chamber output current as a function of 
output power level for the Michigan State University 
TRIGA Mark I nuclear reactor 



CHAPTER III 



DESIGN OF THE REACTIVITY MEASURING DEVICE 



3.1 Formulation of the Problem 



Having presented the equations defining reactivity and its 
relationship to the output current of the ion chamber it now re- 
mains to develop a means by which p can be measured. 

Equation (1.1) expresses the relationship between 0^ 

and I where 
s 



i 

s 



(t) 



V 




(i.i) 



The objective is to process this current so as to produce a voltage 
proportional to the reactor stable period. If this can be accomplished 
effectively then 0^ may be related to p as in Equation (1.2) 
s ince 



0 



P 




( 1 . 2 ) 



Because this relationship is slightly nonlinear as the 
graph in Figure 2.2 illustrates, the value of p corresponding 
to a given 0^ may be taken from this graph and used to calibrate 
the reactivity meter. 
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3.2 Possible Methods of Solution 

Given the desired output p and the input i^(t) several 
methods for determining p will now be considered. 

One method makes use of the circuit in Figure 3.1 where a 
d-c voltage proportional to 1 / 0^ is obtained from i^(t) using 
a logarithmic amplifier followed by a differentiator and filter. 




Meter 



Figure 3.1 Block diagram of a circuit which provides a d-c 
voltage proportional to the inverse of the 
reactor stable period 



In a circuit of this type the logarithm of the ion chamber 
current is first taken such that 

v i (t) ■ - A lo *io i s (t) ■ - A log ioV 0p * - A los ioV^ log i0 e (3 - l > 

The voltage v^(t) is a linear ramp function whose slope 

is inversely proportional to the negative inverse of 0 . Dif- 

P 

ferentiation of this ramp voltage yields a d-c voltage where 



V 2 -ItV') = 



-Kd , . . _ At 

(-A log I - — 

dt s 10 0 9 

P 



l°g 10 e) = 



KA 

9_ 



lo8 io e 



(3.2) 



The first two methods presented were suggested by Dr. L. Giacoletto 
of Michigan State University. The third method was developed by 
the author. 



17 



The result is a d-c voltage inversely proportional to . This 

voltage is then passed through a low-pass filter to eliminate any 

low-frequency noise components passed by the differentiator. 

Finally, due to the slight nonlinear relationship between 

Q and p the voltage V (Figure 3.1) is calibrated on a 

p out 

meter via Figure 3.6 to indicate exact values of reactivity. 

A second method for obtaining a d-c voltage proportional 
to 0 involves the use of an integrator and a quarter-square 
multiplier. The circuit is pictured in Figure 3.2. 



i s <t) 




I Meter 



Figure 3.2 Circuit for obtaining a d-c voltage proportional 
to 9 



In this circuit the input current (i (t)) is first con- 

s 



verted to a voltage 



t/e 

v.(t) = i (t)R = I_R e p 
1 s v 0 v 



(3.3) 



using a buffer amplifier whose input resistance is much greater 



than R . 
v 



This voltage is then integrated to obtain v^(t) where 



-i r e t/e 



^ • 5c | V l (t)d,; - - 2 ^" E < e P - » 



(3.4) 
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The quarter -square multiplier is next used to divide out 
the exponential term in V£(t) such that 

t/e 



v.(t) -I R 0 

V . = -^TTT- = (e P 



out V^(t) 



RC 



- 1 ) 



t/e 



-9 p 

RC 



(3.5) 



1 R e 
0 v 



This resultant d-c voltage is proportional to Q and 

P 

can therefore be calibrated on a meter using the graph of Figure 

3.6 to represent the corresponding values of reactivity. 

A third method for finding Q involves the integration 

P 

of a voltage proportional to the change of the input current. 
Figure 3.3 is a simplified diagram of the circuit. 

V (t) 
c 




Figure 3.3 Circuit for obtaining a discrete voltage propor- 
tional to the slope of i g (t) 

The input current is first converted to a voltage v^(t) 
such that 

t/9 



V, (t) = i (t)R = I R e p 
1 s v 0 v 



(3.3) 



This voltage is then fed into a window comparator whose "GO 11 out- 
put (Figure 3.3) acts as a current sink for V ^ v 1 (t) £ V and 

L U 

as a large resistance whenever v n (t) < V or v- (t) > V . The 

l v £ 1 u 

M G0 n condition will allow the current i (t) to flow in the coil 

c 
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of the reed switch and close its contacts. When this occurs the 



d-c voltage V^. is switched to the input of the integrator and 

is integrated for the period of time that v^(t) remains within 

the bounds of and V . This period of time is proportional 

to 0 and is also proportional to the final integrator output 
P 

voltage. The relationship between V _ and 0 can be found 

out p 

by starting with the output voltage integral 

e 



v = — 

out RC 



{*« 



dt 



(3.6) 



The total integration time (0 ) is that period of time that 
^ v^(t) £ and may be calculated by substituting and 



V into the equation for v- (t) where 

"0 -L 



W ’ \ ‘ R vV 



9./9_ 
l P 



(3.7) 



and 



0 / 9 

v.(0 ) = V = R I e u p 
1 u u v 0 



(3.8) 



0^ and 0^ are the times for which v^(t) is equal to and 

respectively. 

Equations (3,7) and (3.8) are then solved for the time 



increment 0^_ where = ® u “ 6^ and 



V 



e„ = 9, - 9 = 9 (£n 



V 
■tn 



t u 0, V 0 R V I () 



(3.9) 



Using Equation (3.9) for the value of 0^_ in the integral 
of Equation (3.6), the voltage ^ out ma y be written as 



V 



_1 9 (-In — 



In 



out RC p R I. R I. 

v v 0 v 0 



(3.10) 



20 



which is a d-c voltage proportional to the stable reactor period. 

As in the first two methods this voltage may be calibrated 
on a meter using the nonlinear graphical relationship of Figure 3.6 
to indicate the corresponding value of reactivity. 

The above method for measurement of reactivity differs 
from the preceeding two in that the output voltage is not a con- 
tinuous representation of but rather a discrete value obtained 

after integrating a d-c voltage over a discrete time period 0 . 

After consideration of the three alternatives the method 
chosen to measure p was the first one described using a 
logarithmic amplifier followed by a differentiator and filter 
(Figure 3.1). 



3.3 The Circuit and Theory of Operation 

With the reactivity measurement solution method chosen, 
the circuit of Figure 3.4 was designed to produce a d-c voltage 
proportional to 0 . 

The transfer function of the circuit is 



(-R 2 c lS ) x 

V out = V l (s) (R 1 C 1 S+1) (R 2 C 2 s+ 1) (R C f s+1) 



(3.11) 



where (s) is the frequency domain expression for v^(t) given 
as 



V 1 (s) 



■ A lo8 10 ( W 



s 



A log 10 e 



(3.12) 



and A is the logarithmic amplifier voltage constant determined 
approximately by 



• + 15 V 
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A* (• 



R y + 



5000 + R 



R, 



3 ) (26 X 1Q~ 3 ) 



.434 



Volts 



(3.13) 



where is a temperature dependant resistor within the logarithmic 

amplifier module. ^ 

1^ is the log amplifier reference current determined 



pr imar i ly by . 

The concept by which the circuit works is as follows: 

Using the ion chamber current as the input, the logarithm of i^(t) 
is first taken. The output voltage (v^(t)) is a linear ramp 
function since 



i (t) i n t/a 

v i (c > ■ - A lo8 io "T7 ■ - A los io ( r e > 

k k 



a i 0 At . 

- A los io r • r lo8 io e 

k p 



v (t) is then differentiated to produce a d-c voltage where 



(3.14) 
16 



Modular Logarithmic Amplifier Model 4116 (Burr-Brown Research 
Corpo , 1971) , p.3 . 

A more complete discription of the logarithmic amplifier and 
its parameters is given in Appendix A. 

1 6 

Since the input current (i^(t)) is an exponentially increasing 

function and since the circuit (Figure 3.4) used to process this 
function has a limit to the absolute magnitude of the input, i (t) 

can only remain "on" for a discrete period of time Q . Therefore, 

in order to be absolutely correct when defining i g (t) the g ate 

function notation should be used where 

G e t (t ' e t/2> " 1 0 s c s \ 

G e (t ' 0 t/2 ) " 0 t < o, t > e t 

9 is the time of the input pulse. 

Using this notation i g (t) becomes 
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V 2 = 



' R 2 C 1 dt 



57 V c > ■ ' r 2 c i JT ( ' a loe io i; 



At . . 

r log io e) 

P 



R 2 C x A 



los io e 



(3.15) 



This d-c voltage is then filtered to eliminate any low- 

frequency noise enhanced by the differentiator. The resultant 

voltage V is inversely proportional to the reactor stable 

out 

period and approximately proportional to the reactivity such that 



R ? C A R_C i A 

V out = -0“ log io e s p "iT" loS io e (3 * 16) 

P 0 

As previously discussed (Section 2.2) the relationship 
between 0^ and p is nonlinear. Calibration of the meter may 
be accomplished by producing known reactor periods and marking 
the meter deflection with the corresponding values of p via the 
graph in Figure 3.6. 

A comprehensive analysis of the circuit design is given in 
Appendices A, B, C and D. 



3.4 Calibration 

The output voltage in the circuit of Figure 3.4 is a d-c 
voltage proportional to the inverse of 0 and approximately 



t/0 r 



V') - V 



G 9 c <C ' 9 t/2 ) 



If the gate function notation is carried with i g (t) two finite 
strength impulses will appear when v^(t) is differentiated in Equa- 



tion (3.15) due to differentiation of the gate function. 

Although these finite strength impulses are of importance (partic- 
ularly since they show up on the output when i g (t) is switched on 

and off) they shall not be carried in the analysis* It should be under- 
stood, however, that i (t) is on for a finite period of time. This 

s 

period of time (0 t ) is discussed further in Appendix B. 
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proportional to p where 



out 



RC A RCA 

l°g e - p log., e 

0 b 10 H b 10 
P 0 



(3.16) 



As previously mentioned the relationship between 0^ and 

p is nonlinear. Figure 3.6 illustrates this relationship for the 

Michigan State University TRIGA Mark I reactor. 

The range of reactor stable periods for which the meter is 

to be calibrated is 10 sec £ 0 £ 40 sec. This range corresponds 

P 

to a range of reactivity between 40 and 19 cents respectively. 

(F igure 3 . 6) 

Calibration of the reactivity meter is accomplished with 
the apparatus in Figure 3.5. 



1 1 — J 


n — 




11 




- — 


- 



Reactor Power 
Level Indicator 




Figure 3.5 Apparatus used to calibrate the reactivity meter 

The reactor is first made critical with the ion chamber 
current equal to Iq. A control rod is then pulled out a discrete 
distance and the ion chamber current becomes 

t/0 

i (t) = I e P 



( 1 . 1 ) 
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The corresponding d-c voltage on the reactivity meter will 



be 




(3.16) 



This voltage is then marked on the meter. 

In order to mark the reactivity corresponding to ^ out 
the graph of Figure 3.6 is used. As this graph indicates the 
reactor stable period must be known. This value is determined by 
using the reactor power level indicator and a timer to find the 
power level doubling time 0^. The period can be calculated from 
the doubling time since 



9 = 1 o 

p £n 2 



(3.17) 



This process is repeated until the range of p is 



adequately covered. 



p. Reactivity (cents) 
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1 2 3 4 5 10 20 30 40 50 100 

Q , Period (seconds) 

P 

Figure 3,6 Graph showing the relationship between reactivity and 

reactor stable period for the Michigan State University 
TRIGA Mark I nuclear reactor (Courtesy of Gulf General 
Atomic) 



CHAPTER IV 



CONSTRUCTION OF THE REACTIVITY METER 

4.1 Circuit Tests, Performance and Specifications 

The circuit of Figure 3.4 was constructed and tested. 
Figure 4.1 shows several photographs of the instrument. 

Testing of the circuit was accomplished by connecting 
the meter input to the ion chamber via a coaxial cable. The re- 
actor was then made critical and a control rod withdrawn to pro- 
duce a reactor period within the desired range. The range of 
the meter was checked for consistancy and variation due to noise. 
These tests revealed the meter deflection to be approximately 77 
percent for a 10 second period and approximately 24 percent for a 
40 second period thus adequately covering the desired range of 

9 . The maximum variation in meter reading due to noise was less 
P 

than +0.4 percent of full scale. 

The minimum reactor power level for which the circuit was 
designed to operate is ten watts. This corresponds to an ion 
chamber current, Iq = 3nA. Values of input current below this 
level will put the logarithmic amplifier in saturation and there- 
fore cause V to register inaccurately on the meter, 

out 

The upper power level was chosen to be 200 watts resulting 
in an ion chamber current of 60 nA. This value is not an absolute 
upper limit and the meter will continue to read accurately up to 
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Figure 4,1 Photographs of the reactivity meter 
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a reactor power level of 1000 watts (i (t) = 280 na) at which time 

s 

the negative temperature coefficient inherent in the reactor will 
cause the meter to lose accuracy. 

The final specifications of the reactivity meter are listed 
in Table 4.1. 



Table 4.1 Specifications of the reactivity meter 



t/9 p 

Input : Exponentially increasing current. %.i (t) = I e ^ 

s u 



0 10 seconds 

p min 

0 40 seconds 

p max 

-9 

1 . 3 X 10 Amperes 

s min 

-9 

I 60 X 10 Amperes 

s max 

Ou tput reactivity (p) in cents 

p 19 cents 

r min 

p 40 cents 



Accuracy (after calibration) +0.5 percent 



Since the values of p on the meter are calibrated from 
a graph, the only error will be that variation caused by noise on 
the output voltage of the filter (Figure 3.4). 



CHAPTER V 



CONCLUSIONS 



The control rod calibration curve in Figure 1.1 is of the 
type for which the reactivity meter was designed to produced. 
Original design specifications called for the measurement of periods 
between 10 and 40 seconds corresponding to reactivities between 
40 and 19 cents. As the graph indicates, these values cover only 
a small portion of the total rod calibration curve. 

In order to cover the entire range of the calibration 
curve the measured values of rod reactivity must be added in incre- 
ments within the range of the meter. For example: Reference to 

the graph of Figure 1.1 shows the reactivity corresponding to a 
rod position of 300 to be approximately 60 cents. If at this 
position the reactor were made critical and the rod pulled out to 
a position of 350, a reactivity of about 40 cents would be recorded 
on the meter. In order to get the total rod worth at this new 
position this value of reactivity would have to be added to the 
previous total reactivity of 60 cents and would therefore yield 
the new total reactivity of approximately 1 dollar. 

It should also be noted that since the total rod worth is 
determined by the sum of incremental reactivities the error in 
measurement of each incremental rod worth becomes very important 
and the total error increases with each sum. 
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The overall results of the tests made on the instrument 
reveal that the reactivity meter works well for the range of input 
currents and reactor periods for which is was designed. 
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DESIGN OF THE LOGARITHMIC AMPLIFIER 



APPENDIX A 



DESIGN OF THE LOGARITHMIC AMPLIFIER 17 



The logarithmic amplifier used in constructing the circuit 
of Figure 3.4 is a Burr-Brown model 4116 modular logarithmic 
amplifier. Its purpose in the circuit is to take the logarithm 
of the input current i g (t) where 

t/e 

i s (t) - I 0 e <> (1.1) 



The resultant output is a linear ramp function given by 



v 1 (t) = 



i (t) I. t/e 

A S a i 0 P 

- A los io A log io e 



I 

A 1 _fi At . 

■ - A log io ir - T s io e 

k P 



(3.14) 



There are two design parameters for this amplifier. These 
are the voltage scale factor A and the reference current 
Figure A.l shows the logarithmic amplifier with exterior connections. 
The electrical specifications are given in Table A.l. 



The design equations, graphs, figures and specifications in this 
section were taken from the pamphlet Modular Logarithmic Amplifier 
Model 4116 (Burr -Brown Research Corp., 1971), pp. 1-5. 
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Table A.l Electrical specifications of the Burr-Brown modular 

logarithmic amplifier model 4116 (current input only) 



Accuracy 


Accuracy, percent of full scale with current 


source input for 0.4nA ^ I < 400 .. 




Input 


Current Source Input 


. . 400 pA to 400 p,A 


Absolute Maximum Current 




Reference Current Range (1^) 


400 nA £ I £ 40 p,A 


Voltage Scale Factor Range (A) 


. . 2/3 V i A s 10 V 


Output 


Current 




Voltage 




Output Impedance at A = 5 




Stability 


Scale Factor Drift (AA/°C) 


+ 0.0005 A/°C 


Reference Current Drift (AI /°C) 

K 

for 0,4 p,A ^ 1^ £ 1 ^A 


— k 


Input Offset Current Drift (AI g /°C) . . .10 


pA at 25°C, 1 pA/°C 


Input Noise - Current Input 1 pA rms 


; , 10 Hz . to 10k Hz . 


Power Supply 


Rated Voltage 


+15 Vdc 


Supply Drain (Quiescent) 


+ 22 mA 
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R 



x 




Figure A.l Burr-Brown modular logarithmic amplifier model 
4116 with external connections 



A design that takes advantage of the full range of out- 
put voltage is best. The values of R and R required to 

x y 

accomplish this may be obtained by first solving Equation (3.14) 
for the full + 10 volt range where 

I 

V- = -A log = 10 Volts (A.l) 

1 max & 10 I. 

k 

I 

V . = -a log r 222 - = -10 Volts (A. 2) 

1 mm 10 I, 

k 

These equations may be added to solve for with the 

result 



I 



k 



% 

(I I . ) 2 

max min 



(A. 3) 



The corresponding voltage scale factor (A) may next be 
found by substitution of Equation (A. 3) into Equation (A.l) for 
1^. This provides an expression for A where 
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A = 



10 



(A. 4) 



log 



10 



max 



min j 



In Chapter IV it was established that the value of I 



mm 



would be 3nA and I would be 60nA. Substitution of these 

max 

values into the expressions for I, and A yields 



I. = 12.96 nA 
k 

A = 15.75 Volts 



(A. 5) 
(A. 6) 



At this point a problem is encountered. Figure A. 2 shows 
the graphical relationships between 1^, A and their design re- 
sistors. Reference to these graphs and to the specifications in 
Table A.l indicate that the values of 1^ and A in Equations 
(A. 5) and (A. 6) are outside the design limits. 

A compromise must be made to reduce the output voltage 
range in order to put 1^ and A within their respective tolerances. 

For the final design a value of R = 11.3 MQ was selected. This 

to x 

corresponds to a reference current of .5 6p,A. 

Using this value for 1^ the best value for A can be 
calculated by letting the output be 10 volts for 1^ m ^ n * Sub- 
stitution into Equation (A.l) yields 



i mi 

8 10 IT 

k 



(A. 7) 



Reference to the graph in Figure A. 2b indicates a value of 
28 kQ for R will make A equal to 4.5 volts. 

y 

The actual output range of the logarithmic amplifier may 
now be calculated where 
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10 k 100 k 1 M 10 M 100 M 



R (ohms) 

y 




0 20 k 40 k 60 k 80 k 

R^ (ohms) 

(b) 

Figure A. 2 Design curves for the Burr-Brown model 4116 
log amplifier, (a) Reference current vs. R 
(b) Voltage constant vs . R^ 
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V, = 10 Volts 

1 max 



V. 



1 min 



A 1 1 

■ - A los io r. 



max 



-4.5 



log 



10 



60x10 

.56X10 



= 4.36 Volts 



(A. 8) 



(A. 9) 



The design parameters of the logarithmic amplifier may 
be summarized where 



A = 4.5 Volts 
R = 28 . k Ohms 

y 

I k = .56 p,A 

R = 11 .3 M Ohms 
x 

Input Range (i g (t)) 3nA to 60nA 

Output Range (v^(t)) 10V to 4.3 6V 
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APPENDIX B 



DETERMINATION OF THE RAMP FREQUENCY COMPONENTS 



The ideal output of the logarithmic amplifier in the circuit 
of Figure 3.4 is a linear ramp function of the form 

I, 



AT 0 At - 

v i< c > ■ - A io sio r ■ r los io 

k p 



(3.14) 

Differentiation of this function will yield a d-c voltage 



where 



V 2 ' - R 2 C l^ v l (t) ‘ 



- r 2 c i 57 ( ‘ A los io I* - f lc> 8 io e) 

k D 



R 2 C 1 A . 

“1“ l ° 8 10 e 



(3.15) 



In order to obtain the true differential as in Equation 
(3.15) it is necessary to differentiate all frequencies of the 
ramp function v^(t) . This therefore requires the knowledge of 
the range of ramp frequency components . 

The highest frequency components occur when the ramp func- 
tion is on for the shortest period of time. This corresponds to 
a reactor stable period of ten seconds. 

The n on M time for this ramp may be found by substituting 
the limits of i g (t) into Equation (1.1) and solving for 0 t where 



I 

max 



= I . e 
min 



e_/e 

t p min 



( 1 . 1 ) 
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40 



or 



-9 -9 V 10 

10 » 10 ■ 3 » 10 e 



0^ = 10-tn 20 a: 30 seconds 



The Fourier series half -range cosine expansion for 



f(x) = x as a recurrent triangular wave is 



. 18 



c , , k 4k. ttx.I 3nx 1 

f(x) = x = J " ~ ( cos k~ + 9 C ° S ~k~ + 25 C ° S 
TT 



5tpc 



+ . . .+ 



(2n-l) 



(2n-l)nx x 

cos ‘ ) 

k 



0 < x < k 
n = 1,2,.. 



(B.l) 



Substituting the ramp function of Equation (3.14) for 
f(x) the Fourier series expansion for v^(t) is 

. *0 At . . , x 0 A log 10 e 

v (t) -A log — log e -A log — - — 

k p k p min 



40 A 



2 

TT 0 



p min 



, i nt , 1 3rrt , 1 5nt 

1OS 10 e l COS 5; + 9 COS 25 COS ~T t 



, 1 7nt , 1 9nt , 1 lint 

+ — cos h — cos 1 — cos 

49 0 t 81 0 t 121 0 t 



+ . . .+ 



(2n-l) ' 



cos 



(2n-l)nt 



0 < t < Q = 30 sec, 
n * 1 , 2 , • * • 



(B.2) 



The radian frequency components for the minimum duration 
ramp and the respective amplitudes of each (P . ) can be determined 
from Equation (B.2). The values are 



18 

Samuel M. Selby, CRC Standard Mathematical Tables (Cleveland, 
Ohio: The Chemical Rubber Co., 1960), p. 147. 
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uo 



40 A log e 

- = ~ = 0.1048 rad/sec, P = —J 1 = 7.9 

1 0 _ 1 2 ^ 



TT 0 



p mm 



3tt o iia a! p 49fcA 1Og 10 e 

uu, = ~ 0.314 rad/sec, P^ = — — 

9 t 



3.2 

9n 9 



= 0.878 



p mm 



c 40 A log e 

u) 5 = — = 0.524 rad/sec, P 5 = — ^ = °- 316 



25tt 0 



p mm 



7jt 

3 t 



in-, = — = 0.734 rad/sec, P„ = 
7 0 n 



40 t A 3o § 10 e 



7 2 

' 49n 0 



= 0.161 



p min 



q_ 40A log e 

oj q = T 2 - = 0.943 rad/sec, P = — ^ = 0.0975 

9 0 t 9 81 tt 0 

p mm 



uu 



= = 1.15 rad/sec, P 



49 A log e 



10 



11 e 



11 121tt 2 9 



0.0653 



p min 

Since the amplitudes of higher frequency components are 



less than 1/144 of the amplitude of uu^ it is reasonable to omit 



them from the analysis. The value of uo 



11 



radian frequency to be differentiated. 



is therefore the highest 
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DESIGN OF THE DIFFERENTIATOR 



Figure C.l is a magnitude plot for the real and imaginary 
components of an ideal differentiator whose transfer function is 



V . (u>) 
out 

V. (u>) 
in 



= j Kuo 



(C.l) 




uu 

(b) 



Figure C.l Transfer function plot for an ideal differentiator, 
(a) Imaginary part, (b) Real part 



A differentiator circuit with the type of transfer function 
in Equation (C.l) has the disadvantage of enhancing noise. 

If the range of frequencies to be differentiated is limited, 
it is possible to reduce unwanted noise by using a compensated 
differentiator whose transfer function is 



Albert S. Jackson, Analog Computation (New York: McGraw-Hill 

Book Co., 1960), p. 147. 
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where 



V ,(u>) 

out 

V .(to) 
m 



iKu) 



(jKu) + 1) (jKw + 1) 



(C.2) 



and 




„ 2 2 

2K (b 



2 2 2 2 2 
(1 - K u) ) + 4K id 



IRufl - K 2 uj 2 ) 

2 2 2 2 2 

(1 - K u, ) + 4K u) 



(C . 3) 



(C.4) 



These real and imaginary components are plotted in Figure 



C.2. 





Figure C.2 Transfer function plot for a compensated 

differentiator, (a) Imaginary part, (b) Real part 



In Appendix B it was determined that the maximum radian 
frequency to be differentiated in the linear ramp function was 

" 1.15 rad/sec. It is therefore desirable to make 1/K greater 
than 1.15 rad/sec. 

. 20 

Figure C.3 is a circuit diagram of the differentiator. 



20 



RCA Linear Integrated Circuits (RCA, 1967), pp . 75-78. 
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The transfer function may be derived as follows: 
De f ine : 



R, 



Z„ = R„ / /C _ = 



f 2 2 jwR 2 C 2 + 1 



Z = R, + -V 
r 1 juoC, 



V . V. and V „ may be written as 
e’ l 01 y 

(Z + Z )//Z. V (Z //Z.) 

y = y f 01 1 , 01 r 1 

e in Z + (Z £ + Z rt ,)//Z . Z c + Z„ . + Z //Z . 

r f Oi i f Oi r i 



V. = 0 

l 



v 01 = " A 0 (u)) (V e - V 



(C.5) 

(C.6) 



(C.7) 



(C.8) 

(C.9) 



Substitution of V and V. into the expression for 

e l 



Vqi yields 



-V. n A ofa)(Z t+ Z n .)Z, 



01 



(Z £ + Z 0 .)Z. + Z r (Z f + Z 0 . + Z.) + A 0 („)Z 1 Z r 



(C.10) 



If R -» a> the exact expression for the transfer function 



can be written as 
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V .Z . - A (uj) Z .Z c 

out _ Oi i 0 W i f 

V. A n (uj)Z.Z + (Z c + Z n .)Z. + Z (Z. + Z + Z.) 

in u i r ' f Oi l r f Oi i 


(C.ll) 



This expression may be simplified if the following assump- 
t ions are made : 

Assume : 



z . » z //z 

1 f r 


(C . 1 2 ) 


Z 0i <<: Z f 


(C .13) 


Z f « A 0 (t«) 


(C .14) 


Z r « A q (<u) 

„ . L . 21 
Equation (C . 9) may then be written 


(C .15) 


v -z* 

out _ f 

V. Z 

in r 


(C . 1 6) 


When the expressions for and Z^ are substituted 

this equation the differentiator transfer function becomes 


into 


V „u t _ -J R 2 C ," 

V in (jR^^ + l)(jR 2 C 2 u, + 1) 

If R^C ^ = R 2^2 then the magnitude plot in Figure C.2 
will correspond to Equation (C.17) with 1/K = 1/R^C^. 


(C.17) 



For a differentiator error of one percent maximum it is 
necessary to establish a corner frequency about ten times the 
maximum frequency to be differentiated. In other words 



DP DP -10 to =10 (0- - 

^1 1 ^2^2 max 11 


(C .18) 


21 

Ibid., pp. 75-78. 
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where uu^ = 1.15 rad/sec and was determined in Appendix B. 

A trade off of accuracy vs. noise is encountered if Equa- 
tion (C.18) is used to design the differentiator. The prupose of 
the low-pass filter in Figure 3.4 is to attenuate the noise passed 
by the differentiator. It was discovered that for l/R^C^ = 10 
= 11.5 rad/sec the filter could not effectively reduce the noise 
without causing an excessive delay in the d-c response due to its 
long time constant. A decision was made to limit the filter time 
constant to 0.6 second. This meant the corner frequency of the 
differentiator had to be reduced until the noise level was tolerable. 
The experimental resultant corner frequency was 1/R^C^ =3.8 rad/sec 
or about three times as great as the highest frequency to be dif- 
ferent iated . 

The magnitude of the error caused by inexact differentiation 
of the maximum ramp frequency will be small since the amplitude of 
this frequency component is only 1/121 of the major ramp frequency. 
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Jackson, Analog Computation , p. 147. 
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DESIGN OF THE FILTER 



Figure D.l is a schematic of the output filter in Figure 

3 . 4 . 

The filter is a simple RC low-pass filter with a unity 
gain output buffer stage to prevent loading. 




Figure D.l Output filter for the reactivity meter 
The transfer function for the filter is 

v „ (to) , 

= i C D 1 

V. n (u>) jR f C f u>+l ^ 

Since the output voltage of the differentiator is 

a d-c voltage with superimposed low-frequency noise it is desirable 
to make the filter time constant long enough to reduce this noise 
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while at the same time have a reasonable response time to the d-c 
voltage . 

Final tests of the reactivity meter indicated that a filter 
time constant of =0.6 seconds reduced the output noise level 

on the meter to less than + 0.4 percent full scale. This value 
for RfCf also gave a d-c response time of about 3 seconds. 
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SPECIFICATIONS OF THE COMPENSATED ION CHAMBER 



The TRIGA Mark I nuclear reactor at Michigan State University 

contains a Westinghouse type WL-8501 compensated ion chamber. Table 

23 

E.l tabulates the specifications of this chamber. 



Table E.l Specifications of the Westinghouse type WL-8501 
Compensated ion chamber 

Maximum Voltage Between Electrodes (dc) 1000 V 

Operating Voltage for Boron Chamber (dc) 580 V 

Operating Voltage for Gamma Chamber (dc) 0 bo -37 V 

13 

Output Resistance, Signal to Case, Minimum 10 Ohms 

Output Capacitance, Signal to Case, Approximately 130 pF 

2 10 2 

Thermal Neutron Flux Range 2.5x10 to 1.5x10 n/cm /sec 

Thermal Neutron Sensitivity 1.5x10 ^ A/nV 

-13 

Gamma Sensitivity, Compensated 2.5x10 A/R/hr 

Since the output resistance of the ion chamber is greater 
13 

than 10 Ohms the chamber may be treated as a current source. 



Westinghouse WL-8501 (New York: Westinghouse Electric Co., 

1 965) , pp . 1 -2 . 
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